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ABSTRACT 

Growth  stresses  in  amorphous  Si02  scales  formed  during  SiC  fiber  oxidation  were 
calculated.  A  numerical  method  using  Deal-Grove  oxidation  kinetics  and  shear-stress  dependent 
SiC>2  viscosity  was  used.  Initial  compressive  stresses  in  SiC>2  of  -25  GPa  from  the  2.2x  oxidation 
volume  expansion  rapidly  relaxes.  At  >1200°,  viscous  flow  of  amorphous  SiC>2  further  relaxes 
stress  to  negligible  levels.  At  700°  -  900°C,  axial  and  hoop  stress  at  the  GPa  level  persist  in  SiC>2 
near  the  SiC-SiCL  interface.  Radial  expansion  of  the  scale  causes  hoop  stress  to  become  tensile, 
and  axial  stresses  are  driven  to  tensile  values  by  the  Poisson  effect.  These  tensile  stresses  can  be 
>1  GPa  for  thick  scales  formed  at  lower  temperatures  on  surfaces  with  high  curvature. 
Approximate  analytical  expressions  for  growth  stress  are  discussed.  Effects  of  viscosity  variation 
as  well  as  other  assumptions  and  limitations  of  the  calculation  method  are  discussed. 


INTRODUCTION 

Constraint  of  the  2.2  x  volume  expansion  during  oxidation  of  SiC  to  Si02  generates  very 
large  growth  stresses.  Microstructural  evidence  for  these  stresses  exists  for  crystalline  scales. 
High  dislocation  densities  in  crystalline  Si02  near  the  SiC-Si02  interface  are  diagnostic  of  high 
shear  stresses  during  growth  of  new  scale.  [1-2]  Axial  cracks  form  in  the  outer  scale  from  tensile 
hoop  growth  stress. [1,  3]  Tensile  growth  stresses  may  decrease  fiber  strengths,  which  is 
significant  for  structural  applications. 

Growth  stress  has  been  extensively  modeled  for  silicon  oxidation,  which  has  volume 
expansion  similar  to  SiC  oxidation.  [4-5]  Recent  models  recognize  that  flow  at  high  stress  is 
viscoelastic  and  use  the  Eyring  model  for  shear-stress  dependent  viscosity. [6-8]  Radial 
compressive  and  tensile  hoop  growth  stresses  are  predicted  for  oxidation  of  silicon  fibers.  [9- 10] 
Axial  stresses  are  generally  ignored.  [6,  9,  1 1]  For  structural  fibers  this  is  an  important  omission, 
because  axial  residual  stress  affects  fiber  strength.  Tensile  stress  increases  silicon  oxidation 
rates, [12-13]  and  recently  this  has  also  been  demonstrated  for  SiC. [14] 

A  method  to  calculate  all  the  growth  stress  components  anywhere  in  a  SiC>2  scale  formed  by 
oxidation  of  SiC  fibers  is  presented.  The  method  involves  scale  discretization  and  separate 
calculation  for  each  increment.  Calculations  are  done  at  temperatures  from  700°  to  1300°C  for 
amorphous  scales.  Growth  stresses  for  SiC  fibers  with  different  radii  are  examined.  Assumptions 
and  limitations  of  the  method  are  discussed.  Complementary  results  and  a  more  thorough 
description  of  the  background  information  and  method  are  given  elsewhere.  [15] 


METHOD 

A  schematic  of  SiC  fiber  oxidation  with  discretization  of  SiC>2  scale  is  in  figure  1 .  There  are 
two  sources  of  stress.  The  first  is  the  elastic  constraint  of  the  2.2  x  volume  expansion 
accompanying  oxidation.  The  second  is  the  circumferential  expansion  of  old  scale  as  it  is  radially 
displaced  outward  by  formation  of  new  scale  (Fig.  1).  This  creates  tensile  hoop  stress  (ctq)  in  the 
outer  scale,[16]  and  tensile  axial  stress  (ctz)  by  the  Poisson  effect. 
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Figure  1.  Schematic  of  SiC  fiber  oxidation. 

The  large  volume  expansion  raises  concerns  for  growth  stress  modeling.  Elastic  constraint 
of  the  2.2x  volume  expansion  causes  ~25  GPa  compressive  stress  in  SiC>2,  much  larger  than 
stresses  for  which  linear  elasticity  is  valid.  [4]  However,  high  shear  stresses  relax  very  rapidly  to 
values  appropriate  for  linear  elasticity  if  SiC>2  viscosity  is  shear-stress  dependent, [6,  17-18]  and 
discretization  can  confine  incremental  stress  differences  to  values  appropriate  for  linear 
theory.  [7,  19-20]  General  models  of  oxidation  growth  stress  are  coupled  diffusion-reaction  and 
fluid-mechanical  problems. [19]  The  oxygen  diffusion  rate  drives  the  rate  that  stress  in  the  oxide 
is  generated.  Stress,  in  turn,  affects  the  oxygen  diffusion  rates.  The  experimentally  measured 
oxidation  rates  for  SiC  inherently  include  the  effect  of  growth  stress.  The  effect  of  stress  on 
oxidation  rate  will  therefore  not  initially  be  considered  when  modeling  growth  stress. 

The  following  assumptions  are  used  in  modeling  the  growth  stress:  1.  Oxidation  volume 
expansion  is  dilational.  2.  Stresses  resulting  from  constraint  of  oxidation  expansion  are  relaxed 
by  flow  of  Si02  with  shear  stress-dependent  viscosity. 3.  Discretization  of  oxidation  to  small 
increments  allows  use  of  linear  elasticity.  4.  Growth  stress  effects  on  oxidation  kinetics  are  not 
considered.  5.  Stress  relaxation  in  the  SiC  fiber  is  negligible. 

Fiber  Oxidation  Kinetics 


Fiber  oxidation  kinetics  will  not  deviate  significantly  from  flat-plate  geometry  kinetics  until 
the  oxidation  product  for  12  pm  diameter  fibers  is  several  microns  thick.[l,  9,  21]  SiC>2  thickness 
(w)  (Fig.  1)  therefore  obeys  Deal-Grove  kinetics  for  flat-plate  geometry:  [1-2,  22] 


w  =  —A 
2  „ 


4Bt 

1  +-y-  1 

A2 


(1) 


The  parabolic  and  linear  rate  constants  are  B  and  B/A,  respectively.  B  and  A  obey  the  usual 
Arrhenius  relationships.  The  SiC  radius  (b)  after  oxidation  is  (Fig.  1): 


b  =  ^jw2(02  —  O)  +  —  Ow  (2) 

where  b0  is  the  original  fiber  radius  and  £1  is  the  SiC/SiCh  molar  volume  ratio. 

Elastic  Growth  Stress 

The  elastic  stresses  and  strains  are  determined  by  modification  of  a  method  used  for 
sequentially  deposited  coatings. [23]  For  discretized  fiber  oxidation,  the  “deposition”  sequence  is 
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reversed;  the  last  oxide  increment  forms  at  the  SiC-Si02  interface,  and  the  oldest  is  at  the 
surface.  The  system  is  divided  into  the  unoxidized  SiC  fiber,  the  SiC>2  added  from  time  t(i  -1)  to 
t(i)  at  the  SiC-SiCh  interface,  and  the  outer  Si02  increment  added  from  t(0)  to  t(i-l)  (Fie.  1).  The 
axial,  radial,  and  hoop  stresses  in  the  SiC  fiber  at  time  t(i)  are  qvSlC(i),  CT.-SlC(i) ,  and  goSlT(i) 
respectively,  and  in  the  two  SiC>2  increments  are  <jzSl°2(i),  ars,02(i) ,  (?eSl02(i)  and  <jzSl  2(i- 1), 

<jr  °2(i-l) ,  <jeSl°2(i-l),  respectively.  The  effect  of  the  ith  Si02  increment  on  elastic  growth 
stresses  is  found  using  strain  compatibility  equations,  which  is  outlined  in  detail  in  other 
publications.  [15]  Stresses  in  older  increments  are  updated  with  the  stress  values  calculated  for 
the  most  recent  increment. 

Relaxation  of  Elastic  Stress 


The  relaxation  of  the  elastic  stresses  for  all  increments  (j=l  to  i)  in  time  increment  At=t(i)- 
t(i-l)  are  calculated  next.  Shear  stress  relaxation  obeys  a  Maxwell  viscoelastic  model,  [6,  24]  and 
the  Eyring  model  for  shear-stress  (x)  dependence  of  glass  viscosity  (q)  is  used  for  SiC>2.[6,  10, 
24-26]  The  initial  elastic  shear  stress  at  t(i  -1)  is  x0.  The  relaxation  of  x0  to  a  new  value  (x)  in 
time  increment  At  is: 
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where  G  is  the  Si02  shear  modulus, Vc  is  the  activation  volume  for  plasticity  in  SiC>2.  k  is 
Boltzmann’s  constant,  xc  is  the  critical  shear  stress  above  which  plasticity  is  significant  (typically 
~100  MPa),  and  T|0  is  the  stress-free  SiC>2  viscosity  that  has  an  Arrhenius  temperature 
dependence. [27]  x0  is  determined  from  the  principal  stresses  for  all  the  increments  (j=l  to  i)  by 
the  usual  methods. [28] 

Relaxation  of  ae(j)and  az(j)  is  proportional  to  x(j)/x0(j)  and  to  their  difference  with  ar(j), 
which  is  a  boundary  condition,  being  zero  at  the  SiC>2  surface  and  near  zero  elsewhere.  The 
relaxed  values  of  ae(j)and  <rz(j)  (o9(j)'  and  aa(j)')  are  determined  by  solution  of: 
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Radial  Displacement  and  Hoop  Stress  Generation 

Relaxation  expands  the  SiC>2  scale  radially.  The  individual  radial  displacement  of  the  jth 
increment  (ur)  is: 


u 


■CD  = 


n 


5LO. 


fisic(1  +  ezi02(i)  +Eli02(j) 


+  *■1 


(i)ifoi(i)) 


(5) 


where  strains  in  [5]  are  calculated  from  the  stresses  in  [4],  The  total  radial  displacement  of  theith 
increment  is  the  sum  of  the  displacements  of  younger  increments.  This  adds  hoop  strain  (eeSl  2) 
to  outer  layers  as  they  are  forced  to  a  larger  circumference  (Fig.  1).  The  new  hoop  strain  in  each 
increment  is: 
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(6) 


-Si.02 
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£|i02  (0  +  ^  «r(j) 
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b(j  —  1)  —  b{j) 

b(j) 


Recalculation  of  Elastic  Stress  in  the  Scale  and  SiC  Fiber  after  Radial  Displacement 

The  stresses  in  each  SiC>2  increment  are  recalculated  for  eeSl02(j)'  by  solving  the  three  strain 
compatibility  equations  for  the  three  principal  stresses.  Revised  axial  stress  ctzSiC  is  computed 
from  the  force  exerted  by  the  SiC>2  scale,  which  is  the  sum  of  the  axial  stress  in  each  SiC>2 
increment  x  area  of  that  increment: 


°?c(0  = 

i=i 


g|L°2  (0  [2c(i)(w(j)  -  w(j  -  1)  +  w(j)2  -  w(j  +  l)2] 

b(i)2 


(7) 


where  c  is  the  scale  +  fiber  radius.  The  revised  radial  and  hoop  stress  in  the  SiC  fiber  are 
calculated  by  determining  the  net  pressure  (pn)  from  the  sum  of  the  pressures  in  each  annular 
increment: 


i-l 
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(8) 


These  revised  stresses  are  added  to  the  next  increment  i  as  the  program  loops  back  to  equations 
[1  -8].  For  SiC  oxidation  when  w(i)  «  b(i),  the  stress  in  SiC  is  much  smaller  than  that  in  SiC>2. 


RESULTS  AND  DISCUSSION 


Growth  stress  calculations  were  done  using  [1-8]  using  scale  discretization  to  500  layers 
(i=500).  Calculations  were  done  for  oxidation  of  SiC  fiber  for  amorphous  scales  of  500  nm 
thickness  (w)  at  700,  800,  900,  1000,  1100,  1200,  and  1300°C  (Fig.  2).  The  sensitivity  of  these 


Figure  2.  Calculations  of  hoop  stress  (ae  -  blue)  and  axial  stress  (cjz-  red)  throughout  a  500  nm 
thick  SiC>2  scale  on  1,3,  and  10  pm  radii  SiC  fibers,  and  SiC  flat  plates. 
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calculations  to  10x  changes  in  scale  viscosity  (r|0)  and  oxidation  kinetics  parameters  (B0)  are 
shown  in  figure  3.  The  Deal-Grove  oxidation  kinetics  for  SiC  fiber  have  been  reported  for  dry  air 
between  700  and  1300°C,  with  A0  =  6.5  x  10'4  m,  B0  =  1.2  x  10'8  m2/s,  Qa  =111  kJ/mol,  and  Qb 
=  249  kJ/mol.[l-2]  A  full  list  of  all  parameters  used  in  the  calculations  is  published.  [15]  The 
variation  in  growth  stress  with  change  in  fiber  radius  was  examined  by  calculations  using  b0  =  1 
pm,  3  pm,  10  pm  and  “flat-plate”  b0  — >oo  (1  km).  Calculations  for  other  thicknesses  and  for 
crystallized  scales  are  reported  elsewhere.  [15] 


nm 

Figure  3.  Change  in  ae  and  az  growth  stresses  with  change  in  viscosity  (r|0)  or  parabolic 
oxidation  kinetics  (B0)  parameters. 

The  changes  of  cieand  az  throughout  the  scale  with  change  in  b0  at  700°C  to  1300°C  are 
evident  (Fig.  2).  As  fiber  radius  (b0)  increases,  ere  and  az  become  less  tensile  and  more 
compressive.  azis  driven  towards  tensile  values  by  the  Poisson  effect  from  tensile  ae  caused  by 
outward  radial  expansion.  For  a  1  pm  radius  fiber,  tensile  ere  develops  at  the  scale  surface  at  all 
temperatures,  and  is  >3  GPa  at  700°C;  tensile  az  is  >  1  GPa.  For  a  10  pm  radius  fiber,  tensile  ere 
does  not  develop  until  T  >  850°C  and  reaches  a  maximum  of  about  600  MPa  at  1000°C;  tensile 
ctz  does  not  develop  until  ~1100°C  and  reaches  100  MPa  at  that  temperature.  Compressive  elastic 
stress  of  —25  GPa  for  ae  and  az  is  rapidly  relaxed  to  —4  GPa  at  700°C  and  —1  GPa  at  1200°C 
just  1  nm  away  from  the  SiC-SiC>2  interface.  At  T  >1200°C  SiC>2  viscosities  are  so  low  that  stress 
levels  are  insignificant  at  scale  locations  more  than  20  nm  from  the  SiC-SiC>2  interface.  For  flat 
plates  ae=  az,  and  stress  is  always  compressive.  Note  that  a  500  nm  scale  takes  ~108  s  to  form  at 
700°C  and  ~105  s  to  form  at  1000°C,  so  the  latter  is  of  more  practical  interest.  Calculations  for 
crystalline  SiC>2  scale  and  for  several  other  scale  thicknesses  are  reported  elsewhere. [15] 
Crystalline  scale  viscosities  are  much  higher  than  amorphous  scale  viscosities,  so  growth  stresses 
are  also  larger.  The  carbon  introduced  in  SiC>2  scales  during  SiC  oxidation  causes  uncertainty  in 
the  actual  SiC>2  viscosity  (t|0).[15]  Effects  of  a  10x  increase  or  decrease  in  t|0  and  B0  on  growth 
stress  at  1000°C  for  b0  =  3  pm  are  shown  in  figure  3.  A  2x  difference  in  ae  is  evident. 
Differences  between  the  effects  of  r|0  and  B0  on  growth  stress  are  only  significant  for  very  thin 
scales. 

For  the  500  nm  thick  scale  a  “continuous-state”  develops  for  ae  tensile  stress  for  T  >  900°C 
for  b0  =  3  pm  and  T  >  1000°C  for  b0  =  10  pm.  For  thicker  scales  the  continuous- state  region 
develops  at  lower  temperatures.  For  x  <  xc  of  100  MPa  [16]  where  a  stress-free  viscosity  (r|0) 

[17]  is  applicable,  an  analytical  expression  for  the  continuous-state  hoop  stress  [ae(ss)]  can  be 
derived:  [15] 


For  x  >  100  MPa  the  stress  dependence  of  viscosity  is  significant,  and  the  shear  stress-dependent 
viscosity  r\  must  be  substituted  for  r|0.[15] 
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SUMMARY  AND  CONCLUSIONS 


A  method  to  calculate  the  all  growth  stress  components  in  Si02  scales  generated  by  the  2.2 x 
volume  expansion  during  SiC  fiber  oxidation  was  developed.  The  method  can  be  applied  to 
fibers  of  other  materials.  High  compressive  hoop  and  axial  stresses  of  ~25  GPa  are  very  quickly 
relaxed  at  all  temperatures.  Radial  expansion  creates  tensile  hoop  stress  in  the  outer  scale. 

Tensile  hoop  stress  eventually  drives  axial  stress  to  a  tensile  state  by  the  Poisson  effect.  Tensile 
hoop  and  axial  stress  can  reach  values  >  3  and  1  GPa,  respectively  for  oxidation  for  long  times  at 
700°  -  900°C  on  1  pm  radius  SiC  fibers.  At  temperatures  greater  than  1200°C  growth  stresses  are 
quickly  relaxed  to  negligible  levels  by  viscous  flow  of  Si02.  The  accuracy  of  the  growth  stress 
calculation  is  likely  to  be  limited  by  the  accuracy  of  the  amorphous  silica  viscosity.  Tensile 
hoop  stresses  reach  steady-state  values  that  can  be  described  by  analytical  expressions. 
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